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By combining the color dipole model of the nucleon with the assumption that cross sections 
behave asymptotically as In^(s), we are able to describe the data for photon, neutrino and hadron 
interactions with protons at all energies; s is the center-of-mass energy of the interacting particles. 
Specifically, we extrapolate the perturbative QCD calculations into the regime of small fractional 
parton momenta x using a color dipole description of the proton target that guarantees an asymptotic 
In^(s) behavior of all cross sections. The ambiguity of introducing a parametrization associated with 
the dipole approximation is mitigated by the requirement that the saturation of the small-a; structure 
functions produces In^ (s)-behaved asymptotic cross sections, in agreement with the data. The same 
formalism allows us to calculate the cross section for the hadronic pair production of charm particles. 

The results, in particular those for the high-energy neutrino and charm cross sections, are relevant 
for evaluating the sensitivity as well as the background in neutrino telescopes. 


I. INTRODUCTION 

The high-energy behavior of photon, neutrino and pro¬ 
ton cross sections on protons cannot be calculated per- 
turbatively when the fractional momenta x carried by the 
constituents become vanishingly small [1-3]. The struc¬ 
ture functions develop a ln(l/a;) divergent behavior that 
results in a violation of unitarity bounds [4] . It has been 
argued for some time that accelerator and cosmic ray 
data favor a In^(s) behavior of hadronic cross sections 
[5]. In fact, a model-independent analytic extrapolation 
of a In^(s) description of the lower energy data on proton- 
proton total cross sections correctly anticipated [6-8] the 
measurements at the LHC and the Auger cosmic ray ob¬ 
servatory [9-11]. In this paper, we extend this successful 
phenomenological approach to photon and neutrino cross 
sections. 

We present the unified dipole model framework [12- 
17] that describes the behavior of yp, up and pp cross 
sections at high energies and small-a:. Perturbative QCD 
calculations break down at high energy when the proton 
has an increasing number of partons with small fractional 
momenta a;. In a parton picture one can simply think of 
saturation as screening resulting from the fact that the 
increasing number of small-x partons have to be confined 
to a high energy proton of finite size. Asymptotically, the 
proton is a black disk of (mostly) gluons with a radius 
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that increases as Ins. The dipole description of the pro¬ 
ton incorporates color saturation [16]. What is new here 
is that we constrain the parameterization associated with 
the dipole model by the requirement that the cross sec¬ 
tion behaves asymptotically as In^(s) [6, 8]. Our main 
result is that a single dipole parameterization that in¬ 
corporates saturation in this way results in a successful 
description of all high-energy cross section data. 

Our results are relevant to high-energy neutrino de¬ 
tectors [18] whose sensitivity is directly proportional to 
the neutrino cross section. When operating in the PeV 
regime, where IceCube recently discovered a flux of cos¬ 
mic neutrinos, the neutrino cross section can be calcu¬ 
lated perturbatively with an accuracy of better than 5%, 
constrained by measured HERA structure functions [19- 
23]. At EeV energies, relevant for the detection of cos- 
mogenic neutrinos produced in the interactions of cos¬ 
mic rays with background microwave photons, this is no 
longer the case and saturation effects must be included 
in evaluating the sensitivity of ARA [24], ANITA[25], 
ARIANNA [26], JEM-EUSO [27], and LUNASKA [28]. 
Our formalism provides a prediction of the EeV-neutrino 
cross section that, although relying on a (dipole) param¬ 
eterization, is directly supported by a wealth of data. 

The same framework can be used to predict the 
hadronic production of charm particles, for instance by 
cosmic ray-air interactions in the atmosphere [29-33]. 
Above energies of ~ 100 TeV, charmed hadron decay into 
neutrinos is the dominant atmospheric background for 
the detection of cosmic neutrinos. While perturbative 
QCD calculations are less reliable at high energies since 
saturation effects, associated with rndy/s 1, are not 
included, our scheme includes these effects and makes 
definite predictions for the cross section using a dipole 
form that incorporates a wide range of data. We find that 
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FIG. 1: Diagram for the —>■ qqX interaction in the dipole 

picture. 


our charm production cross section is in good agreement 
with previous studies using other dipole model parame- 
terizations [34-36]. The charm production cross section 
and charm energy distribution is of great interest because 
a poor understanding of the charm neutrino background 
interferes with the precise characterization of the cosmic 
flux observed by IceCube [37]. 

The paper is organized as follows: in Section II we 
review the 'y*p cross sections in the dipole formalism; in 
Section III we use a hybrid pQCD and dipole approach to 
calculate the vN cross section; in Section IV we calculate 
the pp ccX and total pp cross sections in the dipole 
formalism. We conclude in Section V. 


II. THE DIPOLE FRAMEWORK FOR 
PHOTON-PROTON CROSS SECTIONS 

In the color dipole picture of the proton [12-17, 38, 39], 
the virtual photon-proton production cross section for 
producing qq pairs, shown in Figure 1, has longitudinal 
(L) and transverse (T) contributions given by 

cr'^ip^gq+xiEj) = J dz (fr\^^^^^^{z,r;Q^)\^adip{x,r), 

where E.y is the energy of the photon and x the fractional 
momentum carried by the quark of the struck proton. 
The quantity ^(z, r; (5^)|^ is the probability that 

7 * fluctuates into a qq pair with transverse separation r 
and longitudinal momentum fraction z and is given by 
[12, 40] 

+ z^yxfier) + mlK^ier)] 

( 2 ) 

(3) 

where Cq is the quark charge, Nc the number of colors, 
= zzQ^ -I- rriq, z = 1 — z, and Kq and Ki the modified 
Bessel functions. 

We have extracted adip from the deep inelastic struc¬ 
ture function ^{XjQ"^) (see, for example, [41, 42]) us¬ 


Parameter 

Value 

Cl 

1.475 X 10"^ 

do 

8.205 X 10-^ 

ai 

-5.148 X 10"^ 

02 

-4.725 X 10"^ 

bo 

2.217 X 10"^ 

bi 

1.244 X 10"^ 

^2 

5.958 X 10-^ 

do 

0.71 


2.82 GeV^ 

zo 

2.4 


TABLE I: Values of our model parameters defined for Eq. 
(5). The do value is obtained by matching Eq. (1), using Eq. 
(5), to the Block et al. [46] structure function, while keeping 
the remaining parameters as fitted by Block et al. [46] from 
HERA data. 


ing the high approximation [43] 

d 

(Xdip{x,r) = — (4) 

with zq = 2.4. Exploiting the convenient parametrization 
of F 2 {x,Q‘^) by Block et al. [44, 45], using their most 
recent result [46] for small x, we obtain: 


cTdipix^r) = do 


7r3|.2(i xY 


aoT 


cizl 


-I- A (oiZg -I- 2B (a2ZQ -I- hif'^ + b2f^B) + 2bof^) 
+ f^B(ai+a2B) + z^A^bi+2b2B) (5) 


-4 - logf ', j ) (6) 

\r^x -I- ZqX) 

B = log (l + I] (7) 

where r = pr and the dipole parameters are given in 
Table I. Our approach is to include primarily the small¬ 
er behavior so an overall normalization factor dq is also 
required [47]. 


Our parametrization of adip is determined at low r by 
the high energy behavior of F 2 that correctly incorpo¬ 
rates saturation. Also, note that the large-r behavior 
does not contribute to high energy cross sections because 
it is suppressed by the wave function in Eq. (1). In fact, 
our procedure makes no attempt at describing F 2 for val¬ 
ues of X > 10“^. In Figure 2 we show the dependence of 
our dipole parametrization on x and r = rQs, which is 
related to the color saturation scale Qs = Qo{xo/x)^^'^ 
[48] . In order to compare dipole cross sections as a func¬ 
tion of T, we fix the saturation scale parameters to Qo = 1 
GeV, xo = 1.642 x IQ-^ and A = 0.2194 as in Ref. [16]. 
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The dipole cross section monotonically increases as x de¬ 
creases for large r. 



At this point we can compute the real photon-proton 
charm production cross section. In this case « 0, 
« 0 and X « {2mq)^f {2mpEj), the threshold value 
that dominates the cross section. We fix niq = rric = 
1.275 GeV, the Particle Data Group value [51]. The re¬ 
sult of this calculation is shown in Figure 4 where it is 
compared to the results of other dipole models. The dif¬ 
ference between our result (similar to using the Soyez 
dipole [40]) and the Goncalves and Machado [35] result 
(using the Munier dipole [49]) comes from a combination 
of the choice of charm quark mass and the form of the 
dipole cross section. The Soyez and Munier dipoles use 
the same lancu-Itakura-Munier form [17], but with dif¬ 
ferent parameters. For example, the values for xq differ 
by two orders of magnitude. 

III. NEUTRINO-PROTON CROSS SECTION 


FIG. 2: Our dipole model for different values of x is shown 
from (dark purple) to 10“® (red) with intermediate lines 

every power of ten. The dipole cross section monotonically 
increases as x decreases for large r. 


In Figure 3 we compare our dipole cross section for 
X = 10“® with other parameterizations in the literature 
]16, 40, 49]; details can be found in the appendix. The 
Golec-Biernat-Wusthoff (GBW) dipole [16] only depends 
on T, while the Soyez dipole [40] also depends ln(l/a;). 
Our model for adip displays approximate scaling behavior 
for large r values where different dipole cross sections 
vary by less than a factor of two for x in the range of 
10“^ — 10“®; see also Figure 2. 
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This work 

GBW (1999) 

- - Munier (2004) 

. Soyez (2007) 



r = rQ, = rQo (f^)*^^ 


FIG. 3: Our dipole model compared with other models in the 
literature [16, 40, 49, 50] at a; = 10“®. Note that at small r 
all the models become color transparent and that at large r 
they only have linear scale differences. 


The total neutrino charged (neutral) current cross sec¬ 
tion for an incident neutrino with energy E^, is given by 



with Qmin = 1 GeV and 


dQ^dx 47ra: \Mf + Q'^ J 

X [Y+FJf+ 2(1-y)Fl±Y_xFn (9) 

Here s « 2mpE^, Y± = 1 ± (1 — y)^ with y = Q‘^/{xs), 
Mi = Mw{Mz) for charged (neutral) current interac¬ 
tion, and (“ —”) refers to neutrinos (antineutrinos). 
To leading order in pQGD = 0, while F^{x, Q^) and 
Fg (cc, Q^) are functions of the parton distribution func¬ 
tions, e.g., for a neutrino charged (GG) and neutral cur¬ 
rent (NG) interaction with an isoscalar target nucleon N 
[ 22 ] 


GG : 


— x{yi -\- u ~\~ d -\- d -\- 2s H- 2b -h 2c) 
xF^ = x(u — u-hd—d-l-2s-l-2b — 2c) 


( 10 ) 


NG : 


FJf = x 


-(L^ + R‘^ + L‘^ + Rd))iu + U + d + d) 


+ 2 + 5 + s + 6) 


xF^ = X 


2 (q — Rt + q — ~ + d' ~ d) 

( 11 ) 

where = sin(0u)) is the sine of the weak mixing an¬ 
gle and the weak couplings are given by = 1 — 
4/3s^, Ld = -1 + 2/3s^, Ru = -4/3s^ and Rd = 
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FIG. 4: Photon-proton charm production as a function of the photon energy is shown for our dipole (upper solid orange) 
and other dipole models in the literature, namely Munier [49] (dashed blue), Golec-Biernat et al.[16] (dot-dashed green), and 
Soyez [40] (dotted pink) dipole models. A similar calculation by Gonclaves and Machado [35] (lower solid gray) is shown as a 
reference. 



FIG. 5: Diagram of i^p interaction in the dipole picture. In 
the charge current interaction V = whereas in the neutral 
current V = Z. 

2/3s^. The antineutrino structure functions for charged 
and neutral current are obtained by replacing q ^ q and 
Fg —>■ —F^. For high neutrino energy, i.e. small x, the 
Fg term does not contribute. 

We evaluate F'^jj. using the dipole formalism with 

^t/l = f f ^ 

q Jo J 

( 12 ) 

where ^{z,r;Q'^) corresponds to the wave function 
for a vector boson {W or Z), of virtual momenta Q^, to 
fluctuate into a qq pair with fractional longitudinal mo¬ 
mentum z and transverse spatial separation r. They are 
computed from the diagrams in Figure 5; in the massless 


quark limit [52, 53] 

K,(z,r;Q^)P = 
|V^|,,(z,r;Q2)P ^ 

X 

|^f,,(z,r;Q2)P = 



Q'^ [zz] [z^ + z^] Fg (er), (13) 


[zz]^KUer), 

[-f'« + + -Rrf] 

[zz] [z^ + z^] Fg (er), 

-^ + -Ru + -^d] 

TT^ 

Q2 [zzfKlier) . 


(14) 


(15) 

(16) 


Following the Henley et al. [23] prescription, we switch 
from the pQCD parametrization of the structure func¬ 
tions and the high energy dipole formalism by using Eqns. 
(10) and (11) for x < xq and Eq. (12) otherwise. 

Our result for the total charged current cross section 
is shown in Figure 6. We used the CTlOnnLO PDFs [54] 
and Xq = 10“®. At low energies our calculation agrees 
with the pQCD calculation, but it incorporates the sat¬ 
uration effect at ultrahigh energies that reduce the neu¬ 
trino cross section for E^, > 10® GeV. The ultrahigh en¬ 
ergy cross section is not very sensitive to the choice of Xq 
between 10“^ — 10“®. Our calculation with xq = 10“® 
can be directly compared to the results of Block et al. 
[46]. As a reference, we show the Cooper-Sakar et al. 
[22] calculation without saturation effects; see also [55- 
57]. In Table II, we have tabulated the charge current 
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cross section for the three calculations shown in Figure function is related to the photon wave function, given in 
6. Eq. (2) by [60] 


EVGeV] 




ccc] 

o 

3 

to 




Sarkar et al. 

Block et al. 

This work 

1 

X 

10® 

4.8 

X 

10-33 

4.31 

X 

10-33 

4.88 

X 

10-33 

2 

X 

10® 

6.2 

X 

10-33 

5.46 

X 

10-33 

6.15 

X 

10-33 

5 

X 

10® 

8.7 

X 

10-33 

7.25 

X 

10-33 

8.07 

X 

10-33 

1 

X 

10® 

1.1 

X 

10-32 

8.87 

X 

10-33 

9.71 

X 

10-33 

2 

X 

10® 

1.4 

X 

10-32 

1.07 

X 

10-32 

1.15 

X 

10-32 

5 

X 

10® 

1.9 

X 

10-32 

1.36 

X 

10-32 

1.42 

X 

10-32 

1 

X 

lOi® 

2.4 

X 

10-32 

1.61 

X 

10-32 

1.65 

X 

10-32 

2 

X 

10i° 

3.0 

X 

10-32 

1.90 

X 

10-32 

1.90 

X 

10-32 

5 

X 

10i° 

3.9 

X 

10-32 

2.33 

X 

10-32 

2.28 

X 

10-32 

1 

X 

10^^ 

4.8 

X 

10-32 

2.69 

X 

10-32 

2.60 

X 

10-32 

2 

X 

10“ 

5.9 

X 

10-32 

3.10 

X 

10-32 

2.96 

X 

10-32 

5 

X 

10“ 

7.5 

X 

10-32 

3.69 

X 

10-32 

3.49 

X 

10-32 


TABLE II: Neutrino charge current cross section values for 
Sarkar et. al [22] pQCD calculation, Block et al. [46], and 
this work. 


IV. PROTON-PROTON CROSS SECTION 

In the dipole model the proton-proton cross section is 
given by [13, 59] 

p— ln(2m,/v's) 

= 2 / dy xig{xi,ii) 

^0 

^ ^gN^qq-\-x{,^2iQ ) : 




Oic 




( 20 ) 


We Hrst calculate the cross section for producing charm 
particles. While such a calculation is not new, our 
dipole calculational framework is directly constrained by 
a wealth of data, and can even accomodate the total 
cross section as we will see next. Atmospheric neutrino 
measurements mostly sample the very forward, large 
Feynman-a;p-, cross section; we therefore show in Fig¬ 
ure 9 our results for the da{pp —)■ cc)/dxF distribution. 
We find that differences resulting from different dipole 
parametrizations are smaller than the spread associated 
with different PDFs. In Figure 10 we have plotted the 
total proton-proton charm production cross section for 
different dipole models using the CTlOnnLO gluon PDF 
with /r = \/4mf and as = 0.33. The high energy dipole 
model evaluations of cr(pp —> ccX) are lower than the 
central perturbative evaluation of the charm pair cross 
section shown in, e.g. [61]. The perturbative calculation 
has large uncertainties associated with the scale depen¬ 
dence and for Ep = 10^° GeV and scales dependent on 
factors of rric, the uncertainty band drops to ct ^ 10“^® 
cm^, the level of the dipole prediction. 

Finally, we can use Eqns. (17-19) to calculate the 
asymptotic proton-proton total cross section. We sub¬ 
stitute [12] Eq. (19) by [63] 

^999 — 2 ■^'I’) -f 0'(lip(^X2 ^ -\- (Jdip{^2^ I*)] (^1) 


where g{xi,p) is the gluon distribution function at the 
scale /X, and Xi ^2 satisfy X 1 X 2 — {2mq)‘^ls with 




T,99 


(z,r;Q2)|2 = 2(iV,-l)|x/.®_^(z,r;Q2)|2 , (22) 


- ^exp(?/) 

Vs 

- ^exp(-?/); 

VS 

see Figure 7. Here agN^qq+x is the partonic cross sec¬ 
tion, which in the dipole picture is given by 


and replace mg —>■ rrig, an effective gluon mass, in evalu¬ 
ating xi and X 2 - The dipole picture correctly reproduces 
the ln^(s)-dependence of the pp total cross section at high 
energies. We match the prediction of Ref. [6] that at 
ultra-high energies the total cross section 

app{s) ~ (0.2817) \nis/{2ml)) mb , (23) 


Xi 

X2 


agN^qq+xix2;Q'^) = j |Vf,(2,r;Q2)|2 

X agqq{x 2 ,Z,v) (18) 

The partonic cross section is directly related to the dipole 
cross section [13, 59]: 

9 

dgqq{x2, z,r) = - [adip{x2, zt) -I- adip{x2, zv)] 

- ^(Tdtp{,X2:V) (19) 

where the different terms correspond to the superposi¬ 
tion of the diagrams shown in Figure 8. The gluon wave 


in agreement with the LHC [64] and Auger [65] mea¬ 
surements; see Figure 11. We fit for nig and /i, while 
using as = as{p), and yield nig = 0.154 GeV with 
= 1.69 GeV^ (^i^ = 1.6 GeV^) for CTIO NLO[66] 
(HERAPDF1.5NLO [58]). We plot the resulting fits in 
Figure 11, where we also shown, for comparison. Block et 
al. [6] calculation; at the highest energies where the cross 
section is dominated by the In^(s) term good agreement 
is found between our calculation and Block et al. 

Because the asymptotic cross sections are dominated 
by gluons only, our formalism predicts that all cross sec¬ 
tions become equal with increasing energy, independent 
of the quantum numbers of the projectile and target. 
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FIG. 6: Neutrino cross sections calculated using HERAPDF NNLO PDFs [58] and our dipole model (orange). For comparison, 
we show a recent calculation done in pQCD by Cooper-Sakar et al. [22] (black), which does not incorporate saturation effects 
and a similar calculation by Block et al. [46]. The solid lines represent the charge current cross sections, whereas the dashed 
line portray the neutral current cross section. 



FIG. 7: Diagram of pp —> qqX for one dipole interaction 
topology, the remainder are shown in Figure 8 . In this picture 
the gluon from the projectile proton has fractional parton 
momenta xi and it is modeled through PDF, while the target 
proton gluon has fractional parton momenta * 2^1 and it is 
modeled through a dipole interaction. 


V. CONCLUSION 

We have used a theoretically motivated but approxi¬ 
mate form of the dipole cross section that extends pertur¬ 
bative QCD calculations into the small-cc regime. Ambi¬ 
guities associated with the small-x behavior of the struc¬ 
ture functions are mitigated by requiring a saturation 
that reproduces cross sections that behave asymptoti¬ 
cally as In^(s), in agreement with data. Using Eq. (5) to 
relate E 2 (x, Q^) to adip{x, r), we have used the structure 
function parameterization of Ref. [46] that reflects this 
asymptotic behavior. Using this formalism, we have cal¬ 
culated the high-energy yp, vp and pp cross sections. Our 
results agree remarkably well with other dipole model cal¬ 
culations, and we predict the very high-energy behavior 
of the neutrino cross section in the EeV energy range as 
well as the hadronic production of charm pairs and to¬ 
tal pp cross section. These are of interest to neutrino 
telescopes and high energy cosmic ray experiments. 



FIG. 8: Possible g ^ qq interactions with a gluon in the 
dipole picture. 
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FIG. 9: Proton-proton charm production differential cross section as a function of Feynman Xi? at incident proton energy 
of 10®GeV is shown for our (orange) and Munier [49] (blue) dipoles. Furthermore, solid lines are calculations done using 
HERAPDF 1.5 [58] while dashed lines uses CTIO NNLO[54]. 



FIG. 10: Proton-proton charm production deep inelastic cross section as a function of the proton energy is shown for our dipole 
(orange), and other dipole models in the literature, namely Golec-Biernat and Wusthoff [16] (green), Munier (blue) [49] and 
Soyez [40] (pink). In all cases, we have used GTIO NNLO PDFs [54] parton distribution functions. Furthermore, a recent 
calculation by Bhattacharya et al. [61] (gray) is shown as a reference as well as recent experimental results from the ALIGE 
collaboration [62]. 
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Here, we specify the dipole models used for comparison 
in the text. In Figure 3, we plot the different dipoles 
together assuming the values given in their respective 
references. 

































































FIG. 11: Best fit total proton-proton cross section calculated using the dipole approximation with CTIONLO [66] (solid orange) 
and HERAPDF1.5NLO [58] (dashed light blue) PDFs compared to Block et al. [6] (dashed gray). Furthermore, we include 
only the log^{s) term from Block et al. fit (solid gray line), which dominates the high energy cross section and is in agreement 
with our calculation; the associated gray band is the one sigma error given by the Block et al. fit. 


A. The GBW dipole 

This dipole cross section defined by Golec-Biernat and 
Wusthoff in Ref. [16] is 


OBW^ 
^dip \ 


where the saturation scale is given by Qsix) = 
Qo {^) and the parameter values have been fixed 
at Qo = 1 GeV, A = 0.277, xq = 4.1 x 10“® and 
CTo = 29.12 mb. 


where Y = ln(l/a;), r = rQs{x), and the factorization 
scale given by Qs{x) = Qo 


f 


\ a 

ln(l - No) 

1 1 — exp 


j In^ 

(1 — No)^~^^Y^ 


b = — "OT-a 


and 


B. The Munier and Soyez dipoles 

This dipole cross sections of Ref. [49] and Soyez [40] is 
of the lancu-Itakura-Munier form [17], 


leff{x,r) 


7s + 


ln(2/r) 

kXY 


where 


a 


Munier / Soyez 
dip 


(x,r) = aoN{T,Y) 


fVo(T/2)2T'«=//(“’’'),r < 2 
1 — exp[—aln^(5T)],r > 2 


The Soyez dipole parameter values [40] are Nq = 0.7, 
7 ^ = 0.738, A = 0.220, xq = 0.163 x k = 9.94, 
and do = 27.3 mb. The Munier dipole parameters [49] 
shown in Fig. 3 are Nq = 0.7, js = 0.627, A = 0.175, 
xq = 0.19 X 10“®, K = 9.94, and do = 37.5 mb. 


N{t,Y) = 
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